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We determine the range of thicknesses and refractive indices for which omnidirectional reflection 
from quasiperiodic multilayers occurs. By resorting to the notion of area under the transmittance 
curve, we assess in a systematic way the performance of the different quasiperiodic Fibonacci mul- 
tilayers. 

PACS numbers: 



INTRODUCTION 



II. QUASIPERIODIC FIBONACCI 
MULTILAYERS 



Photonic crystals are periodically structured dielec- 
tric media possessing photonic band gaps: ranges of fre- 
quency in which strong reflection occurs for all angles 
of incidence and all polarizations. They then behave as 
omnidirectional reflectors, free of dissipative losses. Since 
the initial predictions of Yablonovitch Q and John Q , 
photonic crystals have been attracting a lot of atten- 
tion and a wide variety of applications have been sug- 
gested Q. 

In the one-dimensional case, a photonic crystal is noth- 
ing more than a periodic dielectric structure. Bragg mir- 
rors consisting of alternating low- and high-index layers 
constitute, perhaps, the archetypical example Q. In par- 
ticular, quarter- wave stacks (at normal incidence) are the 
most thoroughly studied in connection with omnidirec- 
tional reflection (ODR) H B @, @, H . 

The introduction of Fibonacci multilayers by Kohmoto 
and coworkers spurred the interest for both possible 
optical applications and theoretical aspects of light 
transmission in aperiodic media [IH . H3 . Il5l Hfl IT^ . In 
fact, the possibility of obtaining ODR in quasiperiodic 
Fibonacci multilayers has been put forward recently pH 

mils mm. 

Underlying all these efforts a crucial question remains 
concerning whether quasiperiodic Fibonacci multilayers 
would achieve better performance than usual periodic 
ones. To answer such a fundamental question one first 
needs to quantify the idea of ODR performance in a 
unique manner that permits unambiguous comparison 
between different structures. Only quite recently a suit- 
able figure of merit has been introduced: the area under 
the transmittance curve as a function of the incidence 
angle 23]. In this paper, we resort to this concept of 
area to rank in a consistent way the ODR characteristics 
of these quasiperiodic systems. 



A Fibonacci system is based on the recursive relation 
So = {if}, Si = {L} and Sj = for 3 > 2 - 

Here H and L are defined as being two dielectric layers 
with refractive indices (n#, ul) and thicknesses (d#, dx,), 
respectively. The material H has a high refractive index 
while L is of low refractive index. The number of layers 
is given by Fj , where Fj is a Fibonacci number obtained 
from the recursive law Fj = Fj-i +Fj_2, with Fq = F\ = 
1. 

In order to properly compare the optical response of 
these systems we will rely on the transfer-matrix tech- 
nique. The transfer matrix Mj for the Fibonacci system 
Sj can be computed as [ll| 



M = M#, Mi = M L , 

= M J _iM J _2, j>2. 
The transfer matrix for the single layer H is 



(2.1) 



cos (3h Qh sin /3h 



M 



H 



1 

QH 



sin Ph cos Ph 



(2.2) 



and a analogous expression for L. Here (3h = 
(2tt/ X)riHdH cos Oh is the layer phase thickness, Oh be- 
ing the angle of refraction, which is determined by Snell 
law. The wavelength in vacuum of the incident radiation 
is A. The parameter qn can be written for each basic 
polarization (p or s) as 



Qh(p) 



uh cos 
n cos Oh ' 



Qh (s) 



ncosO 
tlh cos Oh ' 



(2.3) 



where we have assumed that the layer is imbedded in a 
medium of refractive index n. Henceforth will denote 
the angle of incidence and, for simplicity, the surrounding 
medium we will supposed to be air (n = 1). 
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Let us consider a TV-period finite structure whose basic 
cell is precisely the Fibonacci multilayer Sj. We denote 
this system as [Sj] N and its overall transfer matrix is 



Mf° = (Mj) N . (2.4) 



N 



When the unit cell is 62, the resulting multilayer [S2] 
is periodic. For j > 2, [Sj] N are quasiperiodic. 

The transmittance is given in terms of the matrix 

as El 



(AT) 



(2.5) 



where IIM 



W112 



M denotes the sum of the squares of the 
matrix elements. 

In the theory of periodic systems it is well established 
that band gaps appear whenever the trace of the basic 
period satisfies Q 



|Tt(Mj 



> 2. 



(2.6) 



This should be worked out for both basic polarizations. 
The trace map is a powerful tool to investigate this con- 
dition, especially when the index j is high |24|. In our 
context, it reads as 

Tr(M J+1 ) = Tr(M,) T^M^) - T^M^). (2.7) 

This simple recurrence relation allows us to compute eas- 
ily the band gaps. We quote here the first nontrivial 
cases, namely, the pure periodic system S2 = {LH} and 
the first quasiperiodic one S3 = {LHL}, respectively: 



I cos Pl cos Ph — ^lh sm(3 L sm(3 H \ > 1, 

I cos(2p L ) cospH - A LH s'm(2f3 L ) sinp H \ > 1. 
The function Klh is 

1 / Ql oh 



(2.8) 



2 \q H 



Ql 



(2.9) 



which is frequency independent but takes different val- 
ues for p and s polarizations. However, one can check 
that, irrespective of the angle of incidence, the following 
relation for both basic polarizations holds: 



< 1. 



(2.10) 



Due to the restriction (|2.1U|) . whenever Eqs. (|2.8|) are 
fulfilled for p polarization, they are always true also for 
s polarization. In consequence, the p-polarization bands 
are more stringent than the corresponding s-polarization 
ones IH. 



III. ASSESSING ODR FROM QUASIPERIODIC 
MULTILAYERS 



We first investigate the range of layer thicknesses for 
which ODR exists; that is, when condition ([2.6)) holds 
true for all the incidence angles. Although for the simple 
periodic system S2 analytic approximations are at hand, 
the general problem seems to be very involved and we 
content ourselves with a numerical exploration. 

For definiteness, we fix the refractive indices to the val- 
ues til = 1.75 and uh = 3.35 at A = 10 jam. In Fig. 1 we 
have plotted the zones of ODR for the basic periods Sj 
(with j = 2,3,4,5) in terms of the adimensional thick- 
nesses n^d^jX and n^d^/A. Note that the use of these 
adimensional variables not only simplifies the presenta- 
tion of the results, but, as dispersion can be neglected, 
the results apply to more general situations. 




• LH < 
- LHL [ 

• LHLLH 

= LHLLHLHL 
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n L d L / X 



FIG. 1: Regions where ODR (for p polarization) occurs for 
the Fibonacci systems S2 = {LH}, S3 = {LHL}, S4 = 
{LHLLH}, and 5 5 = {LHLLHLHL}. We have taken 
riL = 1.75 and nu — 3.35 at A = 10 pm. The inset iden- 
tifies the filled ellipses. The marked points correspond to the 
minimum area for each one of the systems. 

The contours of all these regions are approximately el- 
liptical. For every allowed value of n^d^lX there are two 
values of uhcIh / A. This can be traced back to the explicit 
form of Eqs. (|2.8|) for the band gaps. The ellipses for £2 
are the biggest, which confirms that this simple system 
has the best range of ODR in terms of nd/X variables. 
Note also that the usual Bragg solution with layers of a 
quarter- wavelength thick at normal incidence, namely 



n L d L /X = n H d H /X = 1/4, 



(3.1) 



works for S2, but not for the others. 

It is worth stressing the fact that the regions of ODR 
for 62 and 63 are disjoint. This increases the difficulty 
of comparison between these systems. On the contrary, 
all the quasiperiodic multilayers have a significant region 
of common parameters. In fact, from the system Sq on- 
wards, all the elliptic contours are essentially the same 
as for the £5 . 
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FIG. 2: Area under the transmittance curve, denned in 
Eq. (|3.2f) , as a function of til cIl / A and n# (1h / A for the system 
5*2, with the same data as in Fig. 1. 



FIG. 3: Regions of ODR for the same Fibonacci multilayers 
as in Fig. 1 in the plane (til,tih) of refractive indices. The 
curves show the limit of ODR for each stack with the optimum 
thicknesses marked in Fig. 1. 



These regions of ODR are not enough to f ully quantify 
the performance of the multilayer. In Ref. 23] we have 
proposed that, once the materials and the wavelength are 
fixed, the area under the transmittance as a function of 
the angle of incidence 6 



,tt/2 



(JV) 



(0) dO, 



(3.2) 



is an appropriate figure of merit for the structure: the 
smaller this area, the better the performance as ODR. In 
Fig. 2 we have plotted this area as a function of n^d^jX 
and uhcIh/^ for S2. The area has been computed solely 
for the points fulfilling the ODR condition, so the abrupt 
steps give the boundaries of ODR plotted in Fig. 1. How- 
ever, this function varies significantly in the ODR region. 

In fact, for the present case the minimum of this area 
is reached at the point 

n L d L /\ = 0.34305, n H d H /\ = 0.25416. (3.3) 

While the value of n^d^/A essentially coincide with the 
quarter wavelength solution ((3.1)) . n^d^jX differs more 
than 30 % of that solution. 

In Fig. 1 we have marked the points of minimum area 
for each one of the Fibonacci systems Sj. We see the 
strong difference for the periodic and the quasiperiodic 
systems. In fact, for the latter (Sj with j > 3) we can 
summarize the results saying that the optimum area is 
reached approximately at the values of the parameters 



n L d L /\ = 1/8, n H d H /\ = 1/4. 



(3.4) 



In our view, this is a remarkable result: from the principle 
of minimum area, we have consistently derived optimum 
parameters for ODR, which differ a lot from the usual 
solutions found in the literature. 

For the thicknesses giving minimum area, we have cal- 
culated the region in the (nz,,n#) plane for which ODR 



occurs. In Fig. 3 we have plotted the boundary of such 
a region for the same Fibonacci multilayers as before: 
above such curves we have the ODR region. It is again 
the periodic system the first in fulfilling ODR: the on- 
set of the ODR curve is at n# ~ 2.5, in agreement with 
previous estimations [l£j. 

Of course, the optimum parameters for the system Sj 
do not need to be optimum for [Sj]^. To elucidate this 
question, we have computed numerically the values of 
n^^/A and n#<i#/A for different systems containing up 
to 42 layers and for the same indices as before. In Table 1 
we have summarized the corresponding data. We have in- 
cluded only results for the five first periods N = 1, . . . , 5, 
since from [Sj] 5 onwards, all the thicknesses are fairly 
stable, while the area tends rapidly to 0, as one would 
expect from a band gap. We can conclude that the opti- 
mum parameters do not depend strongly on the number 
of layers. 
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FIG. 4: Logarithm of the area computed for the systems [Sj] N 
as a function of the number of layers. 
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TABLE I: Optimum parameters for the systems [Sj] N for dif- 
ferent basic periods Sj and the number of periods N ranging 
from 1 to 5. 



Basic Bandwidth # Periods riLdL/A nudu/X 


Area 




1 


0.34305 


0.25416 


1.01660 




2 


0.33807 


0.22187 


0.47147 


& 0.217 


3 


0.30817 


0.23429 


0.17894 




4 


0.29821 


0.23926 


0.06294 




5 


0.29572 


0.24422 


0.02171 




1 


0.11978 


0.26906 


1.07725 




2 


0.12917 


0.26657 


0.47955 


Ss 0.233 


3 


0.13319 


0.26409 


0.17693 




4 


0.13587 


0.26409 


0.06092 




5 


0.13722 


0.26409 


0.02055 




1 


0.14795 


0.25912 


0.48632 




2 


0.16538 


0.25912 


0.10423 


S 4 0.195 


3 


0.16002 


0.29389 


0.01647 




4 


0.16136 


0.29389 


0.00251 




5 


0.16270 


0.29389 


0.00038 




1 


0.14929 


0.28396 


0.24841 




2 


0.15063 


0.28396 


0.01365 


S 5 0.198 


3 


0.15331 


0.28396 


0.00073 




4 


0.15197 


0.28644 


0.00004 




5 


0.15331 


0.28644 


0.00002 



systems fit into a straight line. On the other hand, the 
periodic system S2 lies on another straight line, but with 
a better slope. That is, for a given number of layers of 
the system (and under the hypothesis of optimum thick- 
nesses), the system [LH] N offers better performance than 
any other. 

Of course, one may think that the bandwidth of these 
systems is different. Sometimes the bandwidth is defined 
at normal incidence, and then it has been argued that 
quasiperiodic systems offer fundamental advantages 20] . 
If we denote by A s hort and Ai on g the longer- and shorter- 
wavelength edges for given ODR bands (of the basic pe- 
riod), it seems more appropriate to define the ODR band- 
width as 

B - ^long ~ A s hort ^ ^ 

2 (Along + A s hort) 

Note that this is the appropriate definition in our case. 
Obviously, the parameters chosen for the purpose of com- 
parison must be the ones giving minimum area; i. e., 
optimum ODR behavior. In fact, we have numerically 
checked that the parameters giving optimum area offer 
also a good bandwidth. In Table 1 we have indicated 
this parameter, confirming again that with the proper 
definition the periodic system is the best. 



IV. CONCLUDING REMARKS 



From previous results for the case of Bragg mirrors, it 
is reasonable to assume that the transmittance of [Sj] N 
tends to zero exponentially with the number of layers. 
To test such an ansatz, we have plotted the area (in a 
logarithmic scale) for all these systems. The results are 
presented in Fig. 4. We think that a simple glance at 
this figure is enough to decide on the performance of the 
quasiperiodic systems as omnidirectional reflectors. 

It is quite clear that all the quasiperiodic systems, ir- 
respective of the index j, behave essentially in the same 
way as far as ODR is concerned. All the points for these 



In summary, we have exploited the idea of minimum 
area to fully assess in a systematic way the performance 
of omnidirectional reflectors. Although quasiperiodic 
systems has attracted a lot of interest due to their un- 
usual physical properties, Bragg reflectors offer the best 
performance, although not at a quarter- wavelength thick 
at normal incidence. We believe that the best feature of 
our approach is that it provides a very clear thread to 
deal with omnidirectional reflection properties in a sys- 
tematic way. Our method is general and can be applied 
to any spectral region. 
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